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Abstract. We give an alternative criterion for ODE equivalence in identical 
edge homogeneous coupled cell networks. This allows us to give a simple proof 
of Theorem 10.3 of Aquiar and Dias, which characterizes minimal identical 
edge homogeneous coupled cell networks. Using our criterion we give a formula 
for counting homogeneous coupled cell networks up to ODE equivalence. Our 
criterion is purely graph theoretic and makes no explicit use of linear algebra. 



1. Introduction 

Coupled cell networks are used to represent systems of coupled dynamical sys- 
tems schematically. Such systems appear either in various biological systems. Net- 
works of eight coupled cells modeling central pattern generators in quadrupeds can 
be used to recover the primary animal gaits [3l |7] . One of the important conclusions 
of the theory of coupled cell systems is that the network itself imposes constraints 
on the possible behaviors of the system even when we lack detailed knowledge of the 
behavior of the cells within the network. A recent application to head movement 
that illustrates the importance of this is |^ . For further applications see [TUj . 

Mathematically coupled cell networks are a subclass of vertex and edge labeled 
directed multigraphs with loops. Vertices with the same label represent multiple 
copies of the same dynamical system. Edge labels represent the type of coupling. 
A compatibility condition is imposed that requires every vertex with a given label 
to receive the same set of coupling types as inputs. As with any class of graphs 
there is a natural notion of isomorphic coupled cell networks induced by bijections 
between the sets of vertices. This representation of coupled cell networks follows 
[8]. An alternative approach is outlined in j^. 

Following Stewart and Golubitsky one may associate to each coupled cell network 
a class of ordinary differential equations that are compatible with the network 
structure, the class of coupled cell systems associated to a coupled cell network. 
As was pointed out in [8] it is possible for non-isomorphic coupled cell networks to 
have the same class of coupled cell systems. In this case we term the two coupled 
cell networks O.D.E. equivalent. We will give a full description of the coupled cell 
systems associated to a coupled cell network for the simple case of identical edge 
homogeneous coupled cell networks. For the definition in the general case see 13]. 

Aguiar and Dias [l] examine the structure of O.D.E. equivalence classes for 
such coupled cell networks. They find a collection of canonical normal forms - a 
collection of networks whose number of edges is minimal within the equivalence 
class. This they term the minimal subclass. 
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In this paper we consider the simplest type of coupled cell networks, the identical 
edge homogeneous coupled cell networks. These are simply directed multigraphs 
with loops where every vertex has the same indegree. Aldosray and Stewart gave 
an enumeration of these networks [2 counted up to isomorphism. 

Using a simpler method specific to the case of homogeneous networks we recover 
Theorem 10.3 of IT] which characterizes the minimal subclass in this case. Further- 
more, we are able to give a recursive formula for enumerating the minimal systems 
with a given number of vertices and edges. 

2. Coupled Cell Systems, Coupled Cell Networks, and O.D.E. 

Equivalence. 

We will deal exclusively with identical edge homogeneous coupled cell networks, 
hereafter referred to simply as networks. 

Mathematically such a network is a directed multigraph where loops are allowed 
and where every vertex has the same in-degree. If the constant in-degree is r we 
will call the network degree r. A directed multigraph consists of a set of vertices V 
and a multiset of edges E with elements mV y-V . A multiset may be thought of as 
a function i? : y x 1/ — > N; we call this function the edge multiplicity Junction. The 
condition that every vertex has the same in-degree, r, is then v €V X^uev -^("' = 
r. 

Given an n cell degree r network G = (V, i?), a choice of finite dimensional phase 
space P ~ M"*, and a function F : PxP'' ^ P such that F{xi;yi, . . . yr) is invariant 
under all permutations of the variables yi, - ■ ■ ,yd, we may produce a vector field 
on P". The vector field for the variable Xi associated to cell i is 

Xi = F{xi;xA'),- . -XAi)) 

where Ji , . . . , Jr are the source cells for the r arcs that terminate at vertex i. The 
complete system is 

ii = F{xi;x (i),...xai)) 
Xn = i^(a;„;x („),.. 

Jl Jr 

The set of such vector fields is a subset of the vector fields on P". 

A vector field obtained from a coupled cell network G by a choice of phase space 
and function F is referred to as a coupled cell system, or an admissible vector field, 
associated to G. We may consider the class of all admissible vector fields for a given 
network G and phase space P. We will denote this class of vector fields by Xq. 

Definition: Two coupled cell networks Gi and G2 are called O.D.E. equivalent 
if there exists a network Gj isomorphic to G2 such that for all choices of phase 
space P 




More prosaically, an n cell degree ri network Gi and an n cell degree r2 network 
G2 are called O.D.E. equivalent if there exists a network G'2 isomorphic to G2 such 
that 



ENUMERATING ODE EQUIVALENT HOMOGENEOUS NETWORKS 



3 



(1) for all choices of phase space P and function Fi : P x P' ^ — >■ P there exists 
a function F2 : P x P^^ P such that for all vertices i 

Fi{x^;x A,),...x A^) ) = F2ix.r,x i,i),...XM ) 
where , • • ■ , jr'i are the source cells for the ri arcs that terminate at cell 

(i) (i) 

i in network Gi and /c^ , . . . , kr^ are the source cells for the r2 arcs that 
terminate at cell i in network 

(2) for all choices of phase space P and function F2 : P x P''^ — >• P there exists 
a function Fi : P x P"^^ P such that for all vertices i 

Fi{xi;xA^),.--XA^) ) = F2{xi;x . . .X (i) ) 

Jl Jti "1 "-T-a 

where jj*'' , . . . , jr'^ are the source vertices for the ri arcs that terminate at 
vertex i in network Gi and , . . . , k).^ are the source vertices for the r2 
arcs that terminate at vertex i in network G'2- 
If we consider P — M. and linear functions Pi and P2 then we obtain the notion of 

linear equivalence. It is shown in [3] that linear equivalence and O.D.E. equivalence 

are equivalent. 

3. Network Operations that Preserve O.D.E. equivalence 

In this section we introduce two operations that can be performed on a network 
that preserve the O.D.E. equivalence class. Since we are dealing exclusively with 
homogeneous networks these operations are a small part of the network operations 
considered in [T] . Both Lemma [l] and Lemma [2] can be deduced from the more 
general arguments in [T], in particular from Proposition 7.4. For completeness we 
give proofs of both Lemma [T] and Lemma [2] using only what is required for our 
simpler case. That we may consider only these two network operations and not 
more general operations is crucial for the results in Section [sj 

Here we give go two simple operations on networks that preserve the O.D.E. 
equivalence class of the network. 

(1) Adding loops: A single loop is added to all vertices in the network. 

(2) /c-Splitting edges: Each edge in the network is replaced by k identical copies 
of the edge. 

Intuitively, it should be clear that these operations preserve the O.D.E. equivalence 
class of the network; however, a formal proof is surprisingly difhcult if one does not 
use the notion of linear equivalence. 

Lemma 1. If network G' is obtained from network G by either of the two network 
operations above then G and G' are O.D.E. equivalent. 

Proof. Using [4j it is enough to prove that the two networks are equivalent when 
the variables Xi are taken to be in M and the function F is taken to be linear. In 
this case we observe that for a degree r network the function P must take the form 

F{x;yi, . . . ,yr) =ax + b{yi H h?>). 

Consider a degree r network. Adding a loop to every vertex we obtain a degree 
r + 1 network. 

Given a function P,. : M' M defined by 

Fr{x;yi, . . . ,yr) = ax + b{yi H h y,-)- 
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we define a function Fr+i : W^^ — > M by 

Fr+i{x;yi, . . ■,yr+i) ^ {a - b) x + b{yi H h yr+i)- 

Clearly we have Fr+i(x] x,yi, . . . , yr) = Fr{x; yi, . . . , yr) and consequently the lin- 
ear vector fields admissible for the degree r network are a subset of the linear vector 
fields admissible for the r + 1 degree network. We can easily go the other direction. 
Given any function F^+i : — M of the form 

Fr+i{x;yi, . . . ,yr+i) ^ax + b{yi H hj/r+i) 

we may define a function Fr :W^R by 

Fr{x;yi,. ..,yr) = {a + b)x + b{yi H h yr)- 

Again we have Fr{x; j/i, . . . , yr) — Fr+i{x; x,yi, - ■ ■ , yr) and consequently we see 
that the two networks have precisely the same set of admissible linear vector fields. 

Consider a degree r network. Performing the edge splitting operation we obtain 
a degree k x r network. Given any function Fr : M. oi the form 

Fr{x;yi, . . . ,yr) = ax + b{yi H h yr)- 

we may define a function F^xr : 1^'^^'^ by 

Fkxrix;yi: ■ . ■ ^ykxr) ^ax+ ^{yi H hykxr)- 

Clearly we have 

/c— times fc— times 

(1) Fkxr{x; yi, . . . ,yi, . . . ,yr, . ■ . ,yr) = Fr{x]yi, . . . ,yr) 

and consequently the linear vector fields admissible for the degree r network are a 
subset of the linear vector fields admissible for the degree k ■ r network. 
Given any function Fkxr ■ R'^^^ — > K of the form 

Fkxr{x;yi,. ■ . ,ykxr) = ax + b{yi -\ h ytxr) 

we may define a function Fr : M'' — > M by 

Fr{x;yi, ...,yr)^ ax + kb{yi -\ h ^r)- 

Again equation ([T]) holds, and consequently we see that the two networks have 
precisely the same set of admissible linear vector fields. 

In both cases we see that the operation produces a new network with precisely 
the same set of admissible linear vector fields. Thus we have that the operations 
preserve the O.D.E. equivalence class. □ 

The operations create a network with a larger degree. However, when a network 
has the required structure, the inverse of these operations may be applied to produce 
a network with a smaller degree. 

Lemma 2. For any identical edge homogeneous coupled cell network G, there exists 
an O.D.E. equivalent network Gm with the following properties: 

(1) At least one vertex has no loops, and 

(2) The greatest common divisor of the multiplicities of the edges is 1. 

We will refer to Gm as a reduced network associated to G. If G is not a reduced 
network then Gm has a lower degree than G. 
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Proof. Let s denote the minimum number of loops on a vertex in G. Consider the 
new network G" formed by removing exactly s loops from every vertex. Clearly G" 
has a vertex with no loops. Since G may be obtained from G" by adding s loops we 
see that G and G' are O.D.E. equivalent. Let d denote the greatest common divisor 
of the edge multiplicities in G". We may form a new network Gm by dividing all 
the edge multiplicities by d. Since G' had a vertex with no loops so does Gm- The 
greatest common divisor of the edge multiplicities of Gm is 1 by construction. Since 
we may obtain G' from Gm by splitting each edge into d edges we see that G' and 
Gm are O.D.E. equivalent by Lemma [l] Thus G and Gm are O.D.E. equivalent 
and Gm has the required properties. □ 

The use of Gm to denote the reduced network is not accidental. We will now 
show that Gm is indeed the unique minimal network in the O.D.E. equivalence class 
of G. Since any network is O.D.E. equivalent to such a reduced network, it suffices 
to show that two reduced networks that are O.D.E. equivalent are isomorphic. 

Lemma 3. IfGi and G2 are reduced network, s and Gi and G2 are O.D.E. equiv- 
alent, then Gi and G2 are isomorphic. 

Proof. Let G'2 be the network isomorphic to G2 which appears in the definition 
of O.D.E. equivalence. We will show that Gi and G'2 are equal. If we take the 
phase space P for the cells to be M and consider linear functions of the form 
F{x, yi, . . . , Ur) = ax + b(yi + ■ • • + yr), then we see that for any choice of ai, bi 
there must exist a2, 62, and for any choice of 02, 62 there must exist ai, 61, such that 

(aild + biA) X = (a2ld + 62S) x 

where A is the adjacency matrix associated to Gi, B is the adjacency matrix 
associated to G2, and x = {xi, . . . ,XnY £ K". Since this holds for all x £ M" we 
must have 



(2) aiM + biA^ a2ld + b2B. 

This matrix condition can be reduced to a system of linear equations of two types: 

(3) fli + biAu = a2 + fo2-Bii 1 <i <n 

(4) biA.j = b2Bij l<i,j <n,i^j 



Since both A and B have a zero on the diagonal they must have some non-zero off 
diagonal entries in order to have the required row sums. Now by Q we see that bi 
and 62 must have the same sign and that Aij 7^ if and only if Bij 7^ 0. 

Since both A and B have at least one zero entry on the diagonal, either there 
must he an 1 < i < n such that An = Bn = or there must exist i j such that 
Aii — but Bii > and Ajj > but Bjj = 0. If we assume that there exists 
^ < i,j < n with i ^ j such that An — but Bu > and Ajj > but Bjj = 0, 
then we obtain 

(5) ai = 02 + b2Bii 

(6) ai + biAjj = 02 

from which we immediately get —biAjj = b2Bii which contradicts our earlier ob- 
servation that bi and 62 must have the same sign. Thus there exists an 1 < i < n 
such that An = Ba = and we can obtain from (jsj) that ai =02. 
Thus we must have 

biAij = 62S.H 
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for all 1 < i, j < n. Now 62 divides biAij for all Since the greatest common 
divisor of the entries of ^ is 1 we must have 62 divides bi. Similarly hi divides b2Bij 
for all i, j. Since the greatest common divisor of the entries of i? is 1, we must have 
hi divides &2- Since hi and 62 have the same sign, we must have &i =62- 

Finally we are able to conclude that A = i? so Gi is equal to G2 as claimed. □ 

4. Examples 

First we show how Figure 1 and Figure 2 of [1] are related using our network 
operations. 



(1) (2) (3) 

Figure 1 . Transforming Figure 1 to Figure 2 of [T] using network 
operations. Edge labels represent edge multiplicities. 

Referring to our Figure [l] notice that network (1) satisfies our criterion for being 
a minimal network. If we split each edge of network (1) into 3 edges then we obtain 
network (2), which is O.D.E. equivalent to network (1). If we now adjoin 2 loops to 
each vertex of network (2), then we obtain network (3), which is O.D.E. to network 
(2) and hence O.D.E. equivalent to network (1). 

Next we apply the results of the previous section to the connected 3 cell degree 2 
networks examined in [9, . They note that up to permutation there are 38 connected 
3 cell degree 2 networks but that 8 of them are O.D.E. equivalent to the lower degree 
networks. Each of these 8 is obtained from one of the 4 minimal connected 3 cell 
degree 1 networks by either adjoining a loop to every cell or by doubling all the 
edges, see Figure [2j 

5. Enumeration 

We begin by outlining the work of Aldosray and Stewart in enumerating homo- 
geneous coupled cell networks. They use the counting result known as Burnside's 
Lemma to enumerate all identical edge homogeneous coupled cell networks with n 
cells and degree r counted up to isomorphism. 

To be explicit let us take V = {1, . . . ,n}. Let us denote the set of all multigraphs 
on V with constant in-degree r by ftn,r- The group of bijections on V is the 
symmetric group on n elements, denoted S„. Each such bijection induces a map on 
fin,r- Thus we have a group action of S„ on fln,r- Two networks are related by S„ 
if and only if they are isomorphic networks. Since we are counting the networks up 
to isomorphism what we actually want to count is the number of distinct S„ orbits 
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Figure 2. The minimal connected 3 cell degree 1 networks and 
their associated connected 3 cell degree 2 networks 



in fln,r- Burnside's Lemma is a tool for counting the number of orbits of a group 
action, it states 



I Orb: 



E |Fixn„,.(ff)|. 

ges„ 



where Fix^^ ^{g) = {w S ^n,r ■ 9 ' ^ — '-^}- If 9 and h are conjugate elements of 
S„ then |Fixn„^(g)| — |Fixn^^ and consequently we may sum over conjugacy 
classes rather than individual elements of S„. Suppose that Ci ...,Cm are the 
conjugacy classes in S„. Let gi be some representative of the conjugacy class Ci. 
We may write our sum as 



(7) 



Orbo„,„(S„)| - |gj E \C^m^nU9^ 



i=l 



There is a bijection between conjugacy classes of S„ and partitions of the integer 
n. Following |2j we will denote a partition of n 



ai • 1 + a2 ■ 2 



+ a„ • n 



by [l"i2"^ . . . n""]. The multiplicative form of this notation is perhaps unfortunate 
but should not cause confusion. The strength of this notation becomes apparent 
when we agree that if ai = then the i"' term in the expression may be omitted. 
Using this notation the 7 partitions of n = 5 may be expressed as follows: 



5 • 1 
3-1 
2-1 
1-1 



[1^] 

1 • 2 [1^2^] 
1 • 3 [123I] 
1-4 [II4I] 



1 • 1 + 2 • 2 
1-2 + 1-3 
1 • 5 



[li22] 
[2I3I] 
[511 
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The set of all partitions of n will be denoted by n„. An element of S5 can be 
associated to each p e n„ as follows: 

[1^] (1)(2)(3)(4)(5) [1122] (1)(23)(45) 

[l32i] (1)(2)(3)(45) [2I3I] (12)(345) 

[I23I] (1)(2)(345) (12345) 
[II4I] (1)(2345) 

Every permutation in S5 is conjugate to one of the permutations that correspond 
to a partition of 5. 

Every permutation a € S„ may be written as a product of disjoint cycles in a 
fashion that is unique up to the order to the cycles. The lengths of these cycles 
form a partition on n called the cycle type of the permutate a. The permutation 
corresponding to a given cycle type is called the normal form of the cycle type. 
Every permutation is conjugate to the normal form of its cycle type. 

Looking at the formula ([t]) we see that it would be advantageous to know the 
size of the conjugacy class associated to a given partition of n. The size of the 
conjugacy class corresponding to [l"i2"^ • • • n""] is 



^ ^ lai2"2 . . . n""ai!a2! ■■■aj.' 

If we consider the partition determines the pattern of parentheses 

[12223I] (_)(_)( )( )( ) 

then n! is the number of ways of writing 1, . . . , rt in the blanks. Observing that we 
can permute each cycle cyclically, that is 

(123) (231) (312) 

are all the same 3-cycle, we must factor out the l"i2"^ • • • n"" possible ways of 
expressing all the cycles. Finally we observe that we may permute cycles of the 
same length freely, so we must factor out the ai\a2l . . . q;„! possible orderings of the 
cycles. 

The main difficulty in enumerating the orbits of S„ lies in determining the size 
of the fixed point set Fixo^ ridi)- We will give the formula for this here and refer 
the reader to the details in [2]. 

Definition: Given p € n„ and s e {1, . . . , n} we may define 



fe=l 

where h = gcd(s, k). 

Clearly ^s,p{z) is analytic about and hence we may write 

CO 

*s,p(z) = ^0r(s,p)2;''. 
r=l 

Theorem 1 (Theorem 8.3 [2]). Let n,r e N \ {0} . Let Hn.r denote the number of 
n cell degree r networks counted up to isomorphism. H^.r is given by 

pen„ k=i 
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1 


2 


3 


r 

4 


5 


6 


1 


1 


1 


1 


1 


1 


1 


2 


3 


6 


10 


15 


21 


28 


n 3 


7 


44 


180 


590 


1582 


3724 


4 


19 


475 


6915 


63420 


412230 


2080827 


5 


47 


6874 


444722 


14072268 


265076184 


3405665412 


6 


130 


126750 


43242604 


5569677210 


355906501686 


13508534834704 



Table 1. The number of n cell degree r networks counted up to 



isomorphism, H, 



We use this theorem to generate Table [T] 

This count however includes disconnected coupled cell networks. From the per- 
spective of dynamical systems we are interested only in the connected identical edge 
coupled cell networks. A disconnected system can be decomposed into a number of 
connected systems. Thus a disconnected n cell network corresponds to a partition 
of n with a„ — i.e. any partition of n except [n^]. 

If we denote the number of connected n cell degree r networks by Kn^r then we 
may enumerate the number of disconnected coupled cell networks as follows 



n 



where 

i^i'n,.r ~^ oft, " 1 



am 

is the number of ways of choosing networks from the Km,r distinct connected 
m cell networks with replacement and where order does not matter. From this we 
obtain 

Theorem 2 (Theorem 10.1 2J). Let n, r e N \ {0} . Let Kn.r denote the number 
of minimal connected n cell degree r networks. We have Ki^r — -ffi,r — 1 Cind for 
n>2 

n-l 

q;=0 



1 



We use this theorem to generate Table [2j 

Now we will use the work of Section 3 to give a recursive formula for enumerating 
the connected minimal coupled n cell degree r networks. 

Theorem 3. Let Mn^r denote the number of minimal connected n cell degree r 
networks. For n > 2 we have Mn i — Kn.i and 



r-l 
s=l 



For n = 1 note that M„ ^ = 0. 
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2 
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5 
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1 


2 


2 


5 


9 


14 


20 


27 


n 3 


4 


38 


170 


575 


1561 


3696 


4 


9 


416 


6690 


62725 


410438 


2076725 


5 


20 


6209 


436277 


14000798 


264632734 


3403484793 


6 


51 


117020 


42722972 


5554560632 


355631996061 


13505066262007 



Table 2. The number of connected n cell degree r networks 
counted up to isomorphism, Kn^r 



Proof: If a connected n cell degree r network is not minimal then it is O.D.E. 
equivalent to a minimal n cell s degree network where s < n. Given a minimal 
n cell degree s network G the question thus becomes how many non-isomorphic 
n cell degree r networks can be obtained that are O.D.E. equivalent to G. We 
have seen that any network G' O.D.E. equivalent to a minimal network G may be 
obtained from G by a combination of adjoining loops and splitting edges (and an 
isomorphism which we may ignore). Let A be the operation of adjoining a root and 
Tk the operation of fc-splitting the edge. Clearly we have Tk oTi = T^i- There is a 
commutation relation between Tk and T^oA = A^ oT^. Using this commutation 
relation we see that any combination of adjoining loops and edge splitting can be 
reduced to a single fc-splitting for some fc > 1 followed by adjoining some number 
of loops. Given that G has degree s and G' has degree r the possible choices of 
fc are constrained by ks < r. Thus there are [r/sj possible values of fc. We then 
adjoin sufficiently many loops to bring the degree to r. 

The number of connected minimal n cell degree r networks is thus given by 



r-1 

E 

s=l 



with the initial condition that Mn.i = ifn,i for n > 2. 
Using this theorem we generate Table [3j 



□ 
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4 


2 


n 3 


4 


30 
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371 


982 


1973 


4 


9 


398 


6265 


55628 


347704 


1659615 


5 


20 


6169 


430048 


13558332 


250631916 


3138415822 


6 


51 


116918 


42605901 


5511720691 


350077435378 


13149391543076 



Table 3. The number of minimal connected n cell degree r net- 
works counted up to isomorphism, Mn,r- 



It is interesting to note that the number of connected minimal 2 cell degree r 
networks for r > 2 is given by (/)(r) where </) is the Euler totient function. The 
appearance of the Euler Totient is explained by the following network diagram: 
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Figure 3. A connected 2 cell degree r network. Edge labels 
represent edge multiplicities. 

In order for a 2 cell network to be minimal at least one vertex must have no 
loops. Without losing generality we may suppose that vertex 2 has no loops. Thus 
vertex 2 must receive r inputs from vertex 1. If we let k, with k < r, denote the 
number of edges from vertex 2 to vertex 1 then we see that vertex 1 must have r — k 
loops. If this network is to be minimal then the three edge multiplicities, r, k, and 
r — k, must be relatively prime. This occurs if and only if r and k are relatively 
prime. For a fixed r the number of 1 < A; < r for which r and k are relatively prime 
is </'(r). Provided that r > 2 we may exclude k = since then r — k = r and all 
edge multiplicities have divisor r and hence the network is not minimal. 

If r = 1 then there are in fact two minimal 2 cell degree 1 networks. 
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Figure 4. The two minimal 2 cell degree 1 networks. 
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